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Abstract 

Tissue fusion eliminates physical voids in a tissue to form a continuous structure and is central 
to many processes in development and repair. Fusion events in vivo, particularly in embryonic 
development, often involve the purse-string contraction of a pluricellular actomyosin cable at the 
free edge. However in vitro, adhesion of the cells to their substrate favors a closure mechanism 
mediated by lamellipodial protrusions, which has prevented a systematic study of the purse-string 
mechanism. Here, we show that monolayers can cover well-controlled mesoscopic non-adherent 
areas much larger than a cell size by purse-string closure and that active epithelial fluctuations are 
required for this process. We have formulated a simple stochastic model that includes purse-string 
contractility, tissue fluctuations and effective friction to qualitatively and quantitatively account 
for the dynamics of closure. Our data suggest that, in vivo, tissue fusion adapts to the local 
environment by coordinating lamellipodial protrusions and purse-string contractions. 


Introduction 

Tissue fusion is a frequent and important event during which two facing identical tissues meet and 
bridge collectively over a gap before merging into a continuous structure (1). Imperfect tissue fusion 
in embryonic development results in congenital defects for instance, in the palate, the neural tube 
or the heart (1). Epithelial wound healing is another illustration of tissue fusion through which a 
gap in an epithelium closes to restore the integrity of the monolayer (2). 

Model in vitro experiments have been developed using cell monolayers to study the different 
stages of healing from collective cell migration to the final stages of closure. In this context, we (3) 
and others (4, 5) have recently demonstrated that, for cells adhering to their substrate, and despite 
the presence of a contractile peripheral actomyosin cable at the free edge, the final stages of closure 
of wounds larger than a typical cell size result mostly from protrusive lamellipodial activity at the 
border. In that case, the function of the actin cable appears to be primarily to prevent the onset 
of migration fingers led by leader cells (6) at the free edge. Gell crawling has also been shown to 
have a major role in tissue fusion in vivo, for example during the closure of epithelial wounds in 
the Drosophila embryo (7). 

However, in physiological developmental situations, there is often no underlying substrate to 
which lamellipodia can adhere to exert traction forces. This is the case, for instance, in neural 
tube formation (8) or in wound healing in the Xenopus oocyte (9). The generally well-accepted 
mechanism in these adhesion-free situations is the so-called purse-string mechanism in which the 
actomyosin cable at the edge of the aperture closes it by contractile activity (10). Note that the 
purse-string and the crawling mechanisms are not mutually exclusive (11) and may be involved 
at different stages of the closure (5, 12). In addition, ’’suspended” cohorts of cells, which do not 
interact with a substrate besides being anchored to a few discrete attachment points, are also 
observed in situations such as collective migration in cancer invasion (13). 

Several experimental studies have documented protrusion-driven collective migration in vitro, 
but the purse-string mechanism has not been thoroughly investigated in model situations. Such 
an analysis imposes to suppress the contribution of the protrusions to closure and, therefore, to 
conduct the experiments on non-adherent substrates. 

In a seminal paper, fibroblast sheets were shown to grow and migrate with their sides anchored 
to thin glass fibers (14). More recent studies extended this observation to keratinocyte monolayers 
or epidermal stem cells bridging between microcontact-printed adhesive tracks (15, 16). However, 
despite recent advances emphasizing the role of tissue remodeling (17), the mechanism of closure 
of a suspended epithelium in the absence of these anchoring sites remains an open question. To 
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Figure 1: Fusion of an epithelium over non-adhering domains. A Schematics of the exper¬ 
iment: the surface (blue) is patterned with non-adhering domains (red) of radius R. Cells (green) 
progressively cover these domains and the cell-free area {A{t)) is dynamically monitored. B 5 x 5 
array of non- adhering domains in phase contrast at the onset of the experiment {R = 42 /im). C 
Time evolution of a single domain {R = 42 /im). 


address this point, we have studied the dynamics of gap closure in an unsupported epithelium in 
which the actomyosin cable and the suspended tissue could not adhere to the substrate. Purse¬ 
string contractility in the absence of protrusions was therefore studied on well-defined mesoscopic 
non-adherent patches within an adherent substrate. 


Results 

We studied the bridging of a monolayer over a well-defined non-adherent gap on adherent glass 
substrates patterned with strictly non-adhesive circular regions of radius R between 5 and 75 /im 
(Figure lA). The surface treatment kept its non-adherent properties for up to 3 weeks in biological 
buffers (18, 19). To ensure we obtained reliable statistics, we worked with arrays of tens of non¬ 
adherent identical domains on which we cultured epithelial Madin Darby Canine Kidney (MDCK) 
cells (Figure IB) (20). Notably, MDCK cell sheets have been previously shown to remain functional 
when suspended over large distances in culture medium (21). Fusion processes in neighbouring 


2 











domains remained independent by imposing a space between each of at least 300 /im, a distance 
larger than the velocity correlation length measured independently for the same cellular system 
(22). Tissue fusion was monitored from confluence (t = 0) up to 4 days. 

Immediately after seeding, the cells adhered on the glass and colonies developed by proliferation. 
The expanding monolayer readily covered non- adhesive domains that had a radius of less than 
10 yum. In these cases, the advancing front edge made no arrest, confirming that cells have the 
ability to bridge over non-adherent defects smaller than their own size (23-26) (Supp. Figure 1). 
At the other extreme, for domains with a radius greater than 70 /im, the monolayer covered only 
the glass surface surrounding the patches (Supp. Figure 1). After several days, we observed the 
development of a tridimensional ”rim" at the boundaries of the domains as already reported (19) 
but no further evolution in the subsequent weeks (19, 27). 

Between these two limiting situations, the monolayer initially surrounded the non-adhesive 
domains and then proceeded to cover them until it eventually fused (Figure 1C, Supp. Figure 1, 
Supp. Video 1). Observations using confocal microscopy at the non-adherent surface / monolayer 
interface revealed the absence of vinculin or paxillin, two proteins associated with focal adhesions. 
This confirmed our basic assumption: the cells did not develop adhesions with the treated surface 
during and after closure (Supp. Figure 2). 

We followed the closure process by monitoring the area covered by cells on domains of various 
sizes over a period of several days. A significant fraction of the domains with radii less than 30 /im 
were already closed when the monolayer reached confluence. Therefore, the mechanism by which 
cells cover these small domains may be different (for example by direct bridging) from the one 
relevant to larger domains. As a consequence, we limited our study to 30 /im < R < bO /im. The 
cell-free area A{t) showed only minor distortions to a quasi-circular shape, allowing us to define 
an effective radius r{t) as r{t) = ^/A{t)/7r (Figure lA). 

It is worth comparing the present experiments with the healing of comparable size wounds of 
the same cell line on homogeneous adhesive substrates in which protrusions at the leading edge 
were shown to be the driving force for closure (3, 4). In both cases, the shape of the hole remained 
relatively circular; in particular, no fingering of the leading edge (6) was observed. However, the 
absence of cell-substrate adhesions drastically slowed down the closing dynamics (typically 30 h 
in the present setting vs. 3 h on a homogeneously adherent substrate for R = 3b /im). Moreover, 
in the experiments described here, the closing was very “noisy” in two respects. First, a given 
hole closed in a seemingly erratic succession of large amplitude retractions and expansions of the 
open area (Figure 2A). In some experiments, we observed closure down to 20% of the initial 
radius, which then re-opened to 50% before eventually closing. However, once the closure was 
fully completed, there was no re-opening (and no indication of a different morphology of the cells 
over the non-adhesive patch compared to the adhesive surface (Figure 1C)). Second, comparing 
several closure events for the same non-adhesive patch size, we observed a very large dispersion of 
the closure times (Supplementary Figure 3, Figure 2C-G). For instance, if R = 3b /rm, the average 
closing time was 44 h and the standard deviation was 18 h (A^ = 150). Because of this large 
dispersion, the entire distributions of the closure times (and therefore meaningful average closure 
times) after 83 h could be accessed only for patches with a radius less than 35 /rm. Unfortunately, 
the development of the above- mentioned 3D rim at the border after typically 4 days prevented us 
from accessing the long-time parts of these distributions for larger domain sizes. 

Given this large variability, we chose to reason in terms of the fraction f{RA) of closed holes 
at a time t for a given initial radius R. This fraction / is plotted as a function of R after 4 days 
(Figure 2H). As previously mentioned, all patches with a radius less than 35/rm closed within 83 
h. By contrast, only a small fraction of the experiments performed at > 55 /im closed in this 
time-frame. As a matter of fact, we never observed the closing of patches whose radius was larger 
than 70 /im. The full dynamic evolution of these fractions is plotted as a heat map in Figure 3A 
for 30 /im < R <bf) /im and 0 < t < 83 h. 

Closure is necessarily a collective effect as cells must form a continuous structure that bridges 
over the non-adherent surface. Indeed, by conducting the experiments in low calcium conditions 
that disrupt cadherin-mediated cell-cell adhesions (19), the efficiency of closure was considerably 
reduced (Supp. Figure 4). 
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Figure 2: Gap closure is a stochastic process. A-B Comparison between experimental (A) 
and theoretical (B) trajectories. The two are visually very close. Note the very large fluctuations 
of the radius during closure. C-E The closure times are widely distributed for a given initial 
radius (bars: experimental values, lines theoretical predictions with 7 = 10 /im^h“^ and D = 
1.6 /im^h“^). Error bars are SEMs. N values are 571 (C), 498 (D), 498 (E). (E) As predicted 
by the model, the ratio {tc)/R^ is approximately constant within the accessible dynamical range 
30 < i? < 36 /im, where long time closure events {tc > 84h) are negligible. The red line corresponds 
to the theoretical prediction (Supp. Eq. (40)) = 2 {^Xd) ~ 0.043h/im“^ with 7 = 10 /im^h“^ 

and D = 1.6 /im^h“^. Error bars are standard deviations. G The coefficient of variation of the 
closure time is approximately constant within the same range, with values consistent with the 

theoretical prediction (Supp. Eq. (46)) indicated with a solid line: ~ ( 3 ^+ 7 ) ^ ~ 

7 = 10 /im^h“^ and D = 1.6 /im^h“^. H Early times are well described by a diffusive process. 
Points are the experimental points resulting from the average of more than 400 trajectories, the 
red line is a linear fit yielding a diffusion coefficient of D = 1.56 ± 0.03 /im^h“^ (error bars are 
SEMs). I Eractions of closed wounds at 4 days (black points). The red line is the theoretical 
fraction (7 = 10 /im^h“^ and D = 1.6 /im^h“^). Error bars are SEMs. 
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Figure 3: The fraction of closed wounds is well described by a stochastic model. A 

Experimental and B theoretical fractions of the closed wounds as a function of the initial radius 
R and time t. The theoretical heatmap B is plotted with the values 7 = 10 and D = 

1.6 /im^h“^ that correspond to the smallest error in the fit of the data (C). The two heatmaps 
are visually very close. C Plane of the error landscape (logarithmic scale). The cross denotes the 
optimum, the red contour bounds the 95 % confidence region. D The closure half-time is the time 
at which half of the wounds have closed. The points are the experimental data (3825 domains) and 
the line is the theoretical time for 7 = 10 /rm^h“^ and D = 1.6 /im^h“^. The fraction of closed 
patches at the end of the experiment is less than 0.5 above i? = 45 jam. Error bars are SEMs. 


No stable lamellipodial protrusions similar to those observed on adherent surfaces were evi¬ 
denced in the present experiments. Moreover, confocal imaging confirmed the presence of a pluri- 
cellular actomyosin cable at the edge of the closing open area (Eigure 4). The contractility of this 
cable was tested with two-photon laser photo-ablation experiments and by inhibiting myosin II 
with blebbistatin. When severed, the cable retracted within a few tens of seconds (Eigure 4B, C, 
Supp. Video 2), indicating that it is under mechanical tension. By contrast, when the epithe¬ 
lium bridging over the non-adherent surface was punctured after closure, the hole did not expand 
upon ablation, indicating that no significant tension is stored in the monolayer itself. These small 
wounds then closed rapidly by developing protrusions presumably on the debris left by the ablation. 
Eurthermore, the addition of blebbistatin almost completely inhibited closure (Supp. Eigure 4), 
whereas the same conditions have been shown to slow down but not halt closure on homogeneous 
adherent surfaces (4). 

Our observations confirm that, as the cells do not interact with the surface in our experiments, 
the contractile pluricellular actomyosin cable along the edge must contract and pull the tissue over 
the adhesion-free surface by a purse-string mechanism. By contrast, the tension in the epithelium 
itself is not a factor in this process. 

To describe these experiments, we wrote the force balance equation at the free edge, on a line 
element of the contractile cable of radius r(t). As ingredients of the equation, we considered (Supp. 
Eigure 5): 1/ a force /cable = ~ 7 /^ due to the line tension 7 of the contractile cable, similar to 
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what has been proposed to describe the shape of single cells anchored to the surface via discrete 
points (28, 29); 2/ a friction force /friction = —/dr/dt where ^ is a friction coefficient encapsulating 
the dissipative processes at the cable and between cells (there is no interaction and hence no 
friction at the monolayer/substrate interface); and 3/ a stochastic force /noise needed to model the 
above-described stochastic effects, such as the wide distributions of closing times (Figure 2C-G, 
Supp. Figure 3) or the very noisy trajectories (Figure 2 A). As puncturing the epithelium did not 
result in opening of the wound, we initially did not include epithelial tension in our description 
(see below). After dividing the force balance equation by the friction coefficient the Langevin 
equation describing the evolution of the radius r{t) reads (Supp. Note 1 ) (30): 

% = ( 1 ) 

where 7 = 7 /^ and \/2Dri{t) is a noise term where the diffusion coefficient D quantifies the 
amplitude of radius fluctuations at the margin. Note that the hypothesis of radial force balance 
is supported by independent force measurements on keratinocytes, showing that, close to the free 
edge, the orthoradial component of the traction stress remained small compared to the radial 
component during closure (17). 

For the sake of simplicity, we further assumed that i/ r]{t) is a Gaussian white noise with 
an autocorrelation function {r]{t)r]{t')) = 6{t — t') and that ii/ 7 , ^ and D remained constant 
(independent of r and t). Note that a constant diffusion coefficient D corresponds to fluctuations 
Act of the epithelial tension a about its average value (zero in the present case): D = Acr^/(2^^) 
(see below). 

During the early stages of fusion, equation (1) reduces to simple diffusion, and the initial mean 
square deviation reads (Supp. Note 4): 

((r(0) - r{t))^)t^o = 2Dt (2) 

The experimental data were in good agreement with this theoretical expression, yielding D = 
1.56 ± 0.03/im^h“^ (Figure 21) and confirming a diffusive behavior of the radius at short times. 
Hence, the cable tension 7 does not contribute to the initial statistics that are fully determined by 
the fluctuations. 

In this framework, the fraction of closed wounds f{R,t) obeyed a backward Fokker-Planck 
equation (30) (see Supp. Notes) : 

which could be solved numerically for a given set of parameters and with boundary conditions in 
accord with our experimental observations: We imposed that r = R was a reflecting boundary (a 
hole never opened beyond the area of the non-adhesive domain) and r = 0 was absorbing (there 
was no re-opening after full closure). The closure time was then the time at which r = 0 was first 
attained (first-exit time). 

A least-squares method allowed us to fit the model to the data over the whole map of the 
fraction of closed wounds (31). Varying R and t at given 7 and D, we minimized the mean square 
standardized error between theoretical frequencies and experimental fractions (see Supp. Note 3 
for the definition of the error function and a full description of the fitting procedure). This fit 
yielded the following estimates of the parameters (Figure 3A): 

7 = 10 /im^h“^ and D = 1.6 /im^h“^ 

where the numbers between brackets give the 95 % confidence interval (Figure 3G). Note that the 
diffusion coefficient is consistent with our previous estimate based on the short-time evolution of 
the closure (Figure 21). With these parameters, the simulated and experimental fractions seemed 
very similar as can be seen in Figure 3A,B. More quantitatively, the particular case of the closure 
half-time at which 50 % of the domains have closed as well as the distributions of closure times for 
various radii were indeed well-described by this set of parameters (Figure 3D, Figure 2G-G, Supp. 
Note 3). Finally, trajectories simulated from equation ( 1 ) (32) with the previously determined 
values of 7 and D closely resembled the experimental ones (Figure 2A,B). 

Altogether, we conclude that our stochastic model provides a self- consistent description of 
the closure dynamics. As the confidence intervals for D and 7 exclude 0, the description is also 
minimal in the sense that none of the components can be removed from the description. 
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Discussion 


We have provided evidence that a cell monolayer can develop over non- adherent surfaces even when 
cells at the edge are not anchored to a substrate to pull on it. This property is the transposition 
at the tissue scale of a single cell’s ability to bridge over defects smaller than their own size. 
Experiments conducted in low calcium conditions or in presence of Blebbistatin show that such 
closures are the result of a collective behavior and that impairing the acto-myosin contractility, in 
particular at the purse-string cable, impairs this process. 

The experimental results are well described by a stochastic model that includes the tension 7 of 
the circumferential actomyosin cable, an effective friction ^ and the amplitude D of the fluctuations 
of the radius reflecting the ones of the epithelial tension. Our theory therefore emphasizes the role 
and function of the purse-string contractility under these conditions. Interestingly, this purse-string 
mechanism is secondary to lamellipodial protrusions when the same MDCK cells are migrating on 
surfaces on which they can develop adhesions (3). As MDCK cells develop a peripheral actomyosin 
cable in both these situations, we conclude that the nature of the substrate on which cells migrate 
controls whether this cable has a regularization or a purse-string function. 

From a typical value of the cable tension 7 = 1 — 10 nN (25, 33), the order of magnitude 
of the (one-dimensional) friction coefficient is ^ = 0.1 — lnN/im“^h. We can compare this value 
with the hydrodynamic two-dimensional friction coefficient measured for the same cells in a similar 
setting but on an adherent substrate (3): ^adh = 10“^ nN/im“^h. The characteristic length defined 
as ^/^adh is typically 100 — 1000 //m, consistent with the correlation length characterizing the 
collective migration of MDCK cells on glass (3, 22). Therefore, it is likely that the friction term 
originates mostly from the monolayer adhering to the glass around the non-adhesive domains. 

An important (and intuitively unexpected) conclusion of our study is that fluctuations actively 
contribute to closing. This is particularly apparent at the onset of closing (short times) where 
fluctuations are actually the dominant term in the closure dynamics (equation (2)). Theoretical 
average closure times can be analytically computed as first-exit times and we obtain for the average 
closure time {tc{R)) = R‘^/2{j-\- D) (Supp. Note 3), which shows immediately that, in a statistical 
sense, a non-zero diffusion coefficient accelerates the closure. The model further predicts that the 
standard deviation of the closure time is proportional to and of same order as the average value 
(Supplementary notes 3,4). Quantitatively, within our limited dynamical range, these predictions 
are borne out by data with the fitted parameters determined previously (Figure 2F,G). Last, we 
validate one of our hypotheses of a white noise in the diffusive term of the Langevin equation (1). 
We used the model, and the fitted parameters, to measure the experimental noise, and checked 
that the autocorrelation function of this noise decays rapidly with time, with a correlation time of 
the order of an hour (Supplementary Figure 6 ). Since experiments are performed over days, this 
confirms that the white noise approximation is indeed appropriate. 

We then checked whether our initial assumption of not including the epithelial tension in 
equation ( 1 ), initially based on tissue photoablation experiments, could be further confirmed. A 
non-zero tension a would add a constant term a = a/^ to the drift coefficient on the right-hand 
side of the Langevin equation ( 1 ). As shown in Supp. Note 2 , analyzing our data with these 
three parameters ( 7 , d, D) yielded a very small (negative) epithelial tension whereas the two other 
parameters retained values close to the ones previously determined (Supp. Figure 7A,B). Therefore, 
the epithelial tension can be safely omitted, suggesting that the time-scales involved in the fusion 
process are sufficiently long to allow elastic stresses to relax (34). Not surprisingly, putting 7 = 0 
in this new equation led to a drastic decrease of the fraction of closed wounds at 83 h in good 
agreement with our experimental observations in presence of Blebbistatin (Supp. Figure 4). 

Next, we asked whether including fluctuations of the cable tension would better describe our 
experimental data than our current hypothesis in which these fluctuations arise only from the 
epithelial tension {D independent of r). Assuming for simplicity that fluctuations in the mono- 
layer tension and in the cable tension are not correlated, we expressed the diffusion coefficient as 
^tot(^) = D D^jr^ where D and are proportional to the (constant) amplitudes of the fluc¬ 
tuations of these two tensions. Fitting model to data within a three-dimensional parameter space 
( 7 ,D,D^) allowed us to conclude that fluctuations in the cable tension could be ignored except 
for values of r much smaller than a cell size (Supp. Notes 3 and Supp. Figure 7C,D). Of note, at 
these very small radius values (r < 1 /rm), the fluctuations of the cable become dominant over both 
the fluctuations of the epithelial tension and the deterministic cable tension (Supp. Note 2). This 
regime corresponds to the very late stages of closure, which are, unfortunately, beyond our exper- 
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imental time resolution. Apart from this regime, the noise term of the stochastic model therefore 
originates from fluctuations of the monolayer tension about its zero mean value. These fluctuations 
result from cell-level dynamics occurring over short time scales, such as rearrangements, divisions, 
or cells being pulled out from the adhesive part. 

From the typical value of ^ determined above, we can also estimate the magnitude of the active 
fluctuations of the epithelial tension: = 2D^‘^ = 10“^ — 1 nN^ /im“^ h. This is, to our 

knowledge, the first experimental measurement of this quantity. We can compare this value to the 
amplitude of tension fluctuations resulting from thermal noise: = 2T)thermai^^- As the 

total friction of the cable over its perimeter is 27ri?^, the amplitude of thermal radius fluctuations 
is T^thermai = ksT/(2 'kR^) whcrc ks is Boltzmann’s constant, T ^ 300 K and R ^ 10/rm. Hence, 
A<^thermai ~ 10“^ —10“^ nN^ /im“^ h : the amplitude of the epithelial tension fluctuations is several 
orders of magnitude larger than that of the thermal tension fluctuations, as expected in these active 
cellular systems (35). 

We end by considering tissue fusion over longer periods of time. A prediction of our model is 
that arbitrarily large wounds will always close, given an arbitrarily large duration, as the stochastic 
process defined by equation (1) yields trajectories that eventually always reach zero (36). Prac¬ 
tically, however, because of the change of behavior induced by the 3D rim that develops at the 
adherent/non-adherent boundary, we stopped our analysis at 83 h. Several arrays were further 
analyzed over longer time periods: at 4 days, the evolution of the fractions of closed wounds was 
still well described by our model (Figure 2H). However, the closure process stopped after that 
time as we observed no evolution between 4 and 7 days (Supp. Figure 8), presumably due to rim 
formation at the edge even though the cells were still active and the border of the wounds was 
still fluctuating. This rim is likely to be the reason why large wound do not close in contrast with 
other cell types such as keratinocytes that are able to close gaps of several 100 /rm (17). 

In conclusion, we have shown that the collective migration of an epithelium can switch between 
two modes, depending on the cells’ affinity for their substrate. Whereas on adhesive surfaces, 
the collective migration is mostly driven by protrusions, our work shows that the purse-string 
mechanism is essential on non-adherent surfaces. Importantly, the active fluctuations of the tissue 
are also crucial and accelerate the closure dynamics. A natural future extension of our work will 
be to elaborate a more complete model where these active fluctuations are combined with the 
viscoelastic rheology of the tissue (17), to get a full description of the closure processes in more 
complex situations including, for example, in vivo tissue fusion in embryonic morphogenesis, or the 
collective migration of cancer cells in fibrillar environments (37). Other cell models may be better 
suited to address these important issues where gaps larger than the ones studied in the present 
work are to be closed. 


Methods 

Preparation of the patterns 

Non-adherent domains of varying radii were micropatterned onto cleaned glass cover slips. The 
glass slides were first uniformly treated by a surface treatment of polyacrylamide and polyethyleneg¬ 
lycol (PEG) on which cells do not adhere (18, 19, 38). Domains whose radii were between 5 and 
75 jiim with a 0.5 /rm increment were defined by photolithography directly on the coating in such 
a way that it remained protected by the photoresist (S1813, Microchem, Westborough, MA) at 
the desired location of the non-adhesive domains (19). Using the photoresist as an etching mask, 
the PEG coating was removed in the photoresist-free areas with an air plasma (Harrick plasma 
cleaner), revealing the underlying glass. The resist was then dissolved leading to PEG-coated do¬ 
mains surrounded by a clean glass surface. The surface treatment was stable for weeks in biological 
buffers (19). 

Cell culture 

MDGK cells (39) were cultured in Dulbecco’s modified Eagle’s medium supplemented with 10 % 
EBS (Sigma), 2 mM L-glutamin solution (Gibco) and 1% antibiotic solution (penicillin (10,000 
units/mL), streptomycin (10 mg/mL)). Gells were seeded and maintained at 37°G, 5 %G02 and 
90 % humidity throughout the experiments. We also used MDGK LifeAct cells (3) for ablation 
experiments (these clones were cultured in presence of geneticin (400/ig/mL)). 
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Figure 4: A contractile pluricellular actomyosin cable localizes at the border of the 
monolayer. A Colocalisation of actin and myosin at the free edge of the closing epithelium. The 
white circle is the underlying domain. B Ablation of the cable (ablation point is figured by the 
white triangle) at t=0 and its subsequent retraction. C Dynamics of retraction of cables on several 
wounds. The red line is the average retraction, the pink area is the standard deviation. 


Blebbistatin (Sigma) was used at a concentration of 50 jaM. Experiments were started in the 
absence of the drug. At confluence, a fraction of the supernatant was pumped out, mixed with the 
drug and re-injected into the well. 

Low-calcium medium (calcium-free DMEM, fetal bovine serum 10%, Penstrep 1%, 50 mM 
calcium) was used to reduce cell-cell adhesion. Experiments were started in regular DMEM and 
the buffer was changed to low-calcium DMEM at confluence. 

Microscopy and data analysis 

The bottoms of Petri dishes or 6-well plates were replaced with patterned glass slides. Cells were 
imaged in phase contrast on an Olympus 1X71 inverted microscope equipped with temperature, 
C02 and humidity regulation (LIS), a motorized stage for multipositioning (Prior) and a Retiga 
4000R camera (Qlmaging). Unless otherwise specified, a lOx objective was used and images 
were acquired every 30 min. Displacements and image acquisition were computer-controlled with 
Metamorph (Molecular Devices). 

Eixed fluorescently marked cells were observed under an upright Imager Z2 spinning disk mi- 
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croscope (Zeiss, Oberkochen, Germany) equipped with a CoolSnapHQ2 camera (Photometries, 
Tuscon, AZ) and a 63x water immersion objective. All acquisitions were controlled using Meta- 
Morph software (Molecular Devices, Sunnyvale, CA). 

Images were processed with the Imaged software or with Matlab (MathWorks, Natick, MA) 
routines. Further analysis was occasionally performed on Origin (OriginLab, Northampton, MA). 

Unless otherwise specified, fractions were computed from at least 100 domains for each size 
measured over at least 4 distinct experiments. 

Laser ablations 

Photoablation experiments were performed on an LSM 710 NLO (Zeiss) microscope equipped with 
a two-photon MaiTai laser and a 40x oil immersion objective. The two-photon laser was used at 
85 % power and at a wavelength of 890 nm. 

Immunoflurocescence 

Cells were fixed in 4% PFA, permeabilized in 0.1% Triton X-100 and blocked in 10% FBS in 
PBS. Vinculin labeling was performed with a mouse monoclonal anti-vinculin antibody (Sigma, 
1:500) and Paxillin labeling was performed with a mouse anti-paxillin antibody (Sigma, 1:500) both 
followed by Alexa 488 donkey anti-mouse (Life Technologies, 1:500). Actin was labeled using Alexa 
546 phalloidin (Life Technologies, 1:1000). Myosin was labeled with rabbit anti-phospho myosin 
light chain (Ozyme, Saint Quentin en Yvelines, France, 1:100) followed by Alexa 488 chicken anti¬ 
rabbit (Life Technologies, 1:1000). Hoescht 33342 (Sigma, 1:10000) was used to mark the nuclei. 
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Supplementary Figure 1: Behaviors of the monolayer according to the radius of the non¬ 
adherent domain. A Domains smaller than a cell size are covered rapidly with no arrest of the 
monolayer. B On domains whose size is comparable with a cell size, the monolayer stops before 
covering them rapidly as it progresses. C For still larger domains (30/rm < R < 65/rm), the 
monolayer surrounds the domains and then covers them by purse string as described in the text. 
D For domains larger than 70 /rm, the monolayer surrounds the domains but cannot cover them. 
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Supplementary Figure 2: The cells do not develop adhesions with their substrate. A-B 

and D-E show no significant signal of the adhesion proteins vinculin (A) or paxillin (D) at the 
basal plane. C,F On the xz sections a thin line void of proteins can be seen over the non-adherent 
surface (figured in red while the adhering surface is blue). The white lines in panels B and D are 
the plane of the sections C and F. A-C: R = 41 /rm; D-F: R = 45.5 /rm. 
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Supplementary Figure 3: Variability of the trajectories. Individual trajectories are noisy but 
also define a very broad distribution of closure times. R = 3b /im. 





Control Low Blebbistatin 
Calcium 


Supplementary Figure 4: Importance of the monolayer cohesion and contractility. Low 

calcium conditions that lead to less cohesive monolayers resulted in a lower fraction of closed wounds 
as did the addition of Blebbistatin that almost halted the closure (Vcontroi = 125, Viow calcium = 25, 
^blebbistatin = 50). R = 40 fiiR. Error bars are standard deviations. 
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Supplementary Figure 5: Schematics of the model. Force balance of a line element at the wound 
margin, for a circular wound of radius r. The three lineic force densities depicted are: the tension 
of the actomyosin cable, a viscous friction force, and a fluctuating force. The initial radius of the 
non-adhesive patch is R. (the cell monolayer is represented in green; the non-adherent surface in 
red). 



Supplementary Figure 6: Noise autocorrelation function. We use the Langevin equation (1) to 
measure the noise term r]{t) = + r ) trajectory, with the fitted parameter values 

7 = 10 /im^h“^ and D = 1.6 /im^h“^. The noise autocorrelation function is obtained by ensemble 
averaging: C{t) = {r]{t)r]{t + r)). It decays quickly, with a correlation time of the order of Ih. 
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Supplementary Figure 7: Model robustness. A,B Influence of an epithelial tension. (A) ( 7 , a) 
plane of the error landscape at D = 1.5 /im^h“^ (logarithmic scale). The cross denotes the 
optimum, the red contour bounds the confidence region. (B) Optimal closure frequency map (linear 
scale), computed for 7 = 7 /im^h“^, D = 1.5 /im^h“^ and a = —0.15 /im^h“^. C,D Influence 
of fluctuations of the cable tension. (C) ( 7 ,!)^) plane of the error landscape ai D = 1.5 /im^h“^ 
(logarithmic scale). The cross denotes the optimum, the red contour bounds the confidence region. 
(D) Optimal closure frequency map (linear scale), computed for 7 = 10 /im^h“^, D = 1.5 /im^h“^ 
= 10 /im^h“^. 
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Supplementary Figure 8 : Long time behavior. Fractions of closed wounds at 4 days (black 
points) and 7 days (blue points). Because of the formation of a peripheral rim, there is no evolution 
after 4 days. Error bars are SEMs. 
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Supplementary Notes 


1 A stochastic description 

1.1 A Langevin equation 

Wound closure dynamics on a nonadhesive patch is a noisy process, with, e.g., large fluctuations 
of the closure time tc at a given radius R. This observation calls for a stochastic description. Since 
the wound is approximately invariant by rotation about its center, force balance on a line element 
at the margin may be expressed as a stochastic equation for the wound radius r{t) 

- + Fi{rR) + F2(r,t) = 0 (4) 

where Fi (r, t) and F 2 (r, t) respectively denote the deterministic and the stochastic component of 
the lineic force density, and ^ is an effective friction coefficient subsuming all dissipative processes at 
play. Since ablation experiments show that the circumferential actomyosin cable is under tension, 
we write 

Fi{r,t) = -^ (5) 

r 

where 7 denotes the line tension of the contractile cable (see Supp. Fig. 5). 

A simple stochastic model of the dynamics may be written as 

f{t) = Di (r, i) + y 2 D 2 (r, t) r]{t) (6) 

where Di and D 2 respectively denote the drift and diffusion coefficients, and 77 (t) is a Gaussian 
white noise with an autocorrelation function {r]{t)r]{t')) = S{t — t'). Dividing the parameter 7 by 
the friction coefficient ^ gives a drift coefficient: 

= (7) 

r 


with 7 = 7 /^. 

Unless explicitly mentioned, we study in the following the stochastic dynamics generated by 
the Langevin equation 

r{t) =V^r]{t) ( 8 ) 

r 

with a constant diffusion coefficient D 

D2{rR)=D (9) 

We observe that collective migration always tends to close the cell-free space, which never opens 
up beyond the area of the non-adhesive patch. The boundary condition ai r = R is therefore 
reflecting. Since closed wounds do not re-open, the boundary condition at r = 0 is absorbing. 
In numerical simulations of the Langevin equation, using a finite time step h [32], the boundary 
conditions are implemented as follows: 

(i) absorbing condition at r = 0 : the simulation stops whenever the radius becomes negative 
r{t + h) < 0; 

(a) reflecting condition at r = R: a radius larger than R^ r{t h) > R is replaced by a radius 
smaller than i?, 2R — r{t h) < R. 

1.2 A Fokker-Planck equation 

Eq. ^ is equivalent to an evolution equation for the transition probability distribution function 
p{r, t \ R^O) between a radius R at the initial time 0 and a radius r at time t [30]. It is convenient 
to write the backward Fokker-Planck equation for p{r, t \ i?, 0): 

^p{r,t\ R,0) = Di{r,t) -^p{r, i | i?, 0) + D 2 (r, t) t\R,0) (10) 
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The probability /(i?, t) that a patch with radius R is closed at time t, or closure frequency, reads 

oR 


f{R^ t) = Prob(tc(i?) < t) = 1 — f p{r, t \ R^O) dr 

Jo 


Integrating (10) over r, we obtain the evolution equation for the closure frequency 


( 11 ) 


( 12 ) 


Since wounds are initially open, the initial condition is /(i?i,0) = 0 for an initial radius R = Ri. 
For all time t, we naturally impose /(0,t) = 1 at i? = 0, and the reflecting boundary condition at 
Ri reads = 0 [30]. 

Using this set of initial and boundary conditions, the numerical resolution of Eq. (12) is per¬ 
formed with the function pdepe of Matlab. For a given set of model parameters, Eq. (12) is solved 
on the interval R G [ 0 ,i?i]. The value at i?i, fth{RiR) = f{R = RiR) can then be compared with 
the experimentally measured fraction fexp{RiR) of wounds of initial radius Ri closed at time t. 

From Eq. 0 , /(i?, t) is also the cumulated distribution function of closure times. As a 
consequence, the distribution function P{tc \ R) of the closure times tc of patches of radius R is 
obtained by differentiating / with respect to time: 


P{R\R) = ^{R,R) 


(13) 


2 Parameter fitting and model selection 

2.1 A least-squares method 

In this section, the theoretical closure frequencies ft}^{RR \ 7 , are computed numerically for a 
given set of parameter values ( 7 ,!)). Initial and boundary conditions for the evolution equation 




(14) 


are given in Sec. 1.2 Varying (RR) at fixed ( 7 , we calculate the mean square standardized 
error between theoretical frequencies and experimental fractions [31]: 




]_ t I 7 , D) — /exp(^, t)Y 

N ^ 


R,t 


(R,t) 


(15) 


where N is the total number of data points, and variance of the fraction of patches 

of radius R closed at time t, as measured over Vexp(i?) experiments: 

fexp{RR) (1 - fexp{RR)) 


T .(^-0 = 


A^exp(-R) 


(16) 


The experimental fractions /exp(i?, t) are measured from a total of 6625 patches, with radii ranging 
from 31 to 49.5 fim with a 0.5 fim step, over a total duration T = 83h, and with a time resolution 
At = Ih. In practice, we define Nr = 13 bins of width A/t =1.5 /rm and Nt = 84 time points per 
radius, so that N = Vr Nt = 1092. 

The error landscape is shown in Fig. 3C. The minimum of the mean square error is achieved 
for 


Dr. 


corresponding to a minimal mean square error value 

77i2 _ Z7'2/~ _ Z7'2/~ . 



(17) 


(18) 

7 min 5 D = -Drain) — 7.6 

(19) 


at the bottom of a well-defined single well. The optimal parameter values ( T7p8 ) yield the best 
agreement with experimental data (compare Figs. 3A and 3B). 
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At a given radius R and time t, the 95 % confidence interval of the optimal theoretical frequency 
/th‘" = I 7min,-Dmin) reads - 1.96cr/^^p, + 1.96where is a proxy for 

the standard deviation of Substituting the upper bound of the confidence interval for 

into Eq. yields an upper bound of the mean square error 


El 


^ ^ ^ {f^'-{R,t) +1.96* t) - /e.p(i?, 

N ^ a'i 

R,t /exp ' 


^min + 2 * 1.96 * 


{R,t) 

1 ^ /th"(^/) -/exp(i?/) 

N 2 ^ 


R,t 


0'/exp(-R/) 


+ 1.96^ 


( 20 ) 

( 21 ) 


We define the confidence region for the parameters ( 7 , D) by the domain within the level contour 
= E\ = 15.8, see Fig. 3C. Conservative estimates of confidence intervals on 7 and D (in 
brackets) are finally obtained by inscribing this contour within a rectangle: 


7 = 10 //m^h ^ [6,13] 

D = 1.6//m^h-^ [0.5,3.9] 


( 22 ) 

(23) 


The value of the dimensionless ratio allows to determine whether drift or diffusion dominate 

rCUi 


the dynamics. Since 


D 2 


RDi 


= k=0.16 [0.04,0.65] 

7 


(24) 


we conclude that drift dominates, but that diffusion cannot be neglected. 


2.2 Influence of an epithelial tension 

Since tissues may quite generally be under compression or under tension, we first tested the robust¬ 
ness of our results by taking into account an epithelial tension a of unknown sign. Force balance 
is modified: a reduced tension coefficient a = a contributes to the drift coefficient. Following 
the procedure given in Sec| 2 . 1 [ here based on the Langevin equation 


r{t) = -- + d VWr]{t) 
r 

the optima and confidence intervals for the three unknown parameters ( 7 , d, D) are 

7 = 7//m^h“^ [0,26] 

D = 1.5//m^h-^ [0.4, 3.9] 

d = —0.15/imh“^ [—0.65, 0.65] 


(25) 


(26) 

(27) 

(28) 


for a minimal value of the error E^^^ = d, D) = 7.2 (see Supp. Fig. 7A). 

Strikingly, zero belongs to the confidence interval for d, and the level of agreement between 
predictions and data is not improved when taking into account the (small, negative) optimal value 
(28) (compare Supp. Fig. 7AB with Fig. 3AB). Further, the revised estimates ( [26p7| ) for 7 and D 
are consistent with the previous confidence intervals ( 22p3 ). We conclude that the closure fraction 
data is consistent with the absence of measurable tension in the monolayer (a = 0). Indeed little 
to no retraction was observed when performing laser ablation in the monolayer away from the 
margin. 


2.3 Influence of fluctuations of the cable tension 

The diffusion coefficient D is proportional to the amplitude of fluctuations of the epithelial tension 
about its (zero) average value. We also tested the robustness of our results by taking into account 
possible fluctuations of the the cable tension, of amplitude > 0 . Assuming for simplicity that 
fluctuations in the cable tension and in the epithelial tension are uncorrelated, we express the 
diffusion D 2 {rR) as 


D2{r, t) 




(29) 
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and study the modified Langevin equation 


f{t) 



2 (d+^),{() 


( 30 ) 


(31) 

(32) 

(33) 


interpreted according to Ito’s rule. This model leads to the estimates 

7 = [4,12] 

D = 1.5[0,2.5] 

= 10Mm'‘.h“M0,870] 

for a minimal value of the error = min.y^£)^£)^ E'^{^,D,D^) = 7.5. We emphasize that: (i) 

the optimal values in ( M|32| ) are consistent with the confidence intervals ( |22p3l ); (ii) di non-zero 
value of ( [3^ has little influence on the level of agreement between theoretical and experimental 
closure fractions (compare Supp. Figs. 6C-D and Figs. 3A-B). 

Taking into account fluctuations of the cable tension allows to define two critical radii: 



(34) 


Fluctuations of the cable tension dominate fluctuations of the epithelial tension below R^\ and 



dominate the deterministic cable tension below Using (3l|33), we find R^^ 3/rm and 

^ 1 /im: this suggests that cable tension fluctuations may dominate the very late stage of the 

closure process. 

Conversely, we find that cable tension fluctuations are negligible except near closure, and con¬ 
clude by selecting the most parsimonious model, Eq. <§ , which we will use in Sec. to define and 
compute additional quantifiers of closure dynamics. 


3 Statistical quantifiers of closure dynamics 

3.1 Closure half-time ti/2 

We first consider the closure half-time ti/ 2 , defined as the time needed to close half of the wounds 
for a given initial radius R: 




(35) 


Experimentally, ti /2 becomes larger than the total duration of the experiment ti /2 > 83 h above 
a radius R = 44.5 /im. In Eig. 3D, we compare our measurement of the half-closure time up 
to R = 49.5/im with the outcome of numerical simulations for the optimal parameters ( T7p8 ). 
Experimental error bars are obtained from the maximum max (t+ — ti/ 2 , ti /2 — t“) where the times 


are 


defined by /exp(R,t^) — ^ ^ n ^iR) \ expected by comparing Eigs. 3A-B, where 

the same closure fraction data and theoretical frequencies are plotted as a heat map, experimental 
and theoretical values agree very well in Eig. 3D. 


3.2 Distribution of closure times 4 

Eigs. 2C-E compare experimental and theoretical closure time distributions for radii R = 34, 40 
and 46 /im. Experimental values are computed, within AR =1.5 /im and binned over time intervals 
of At = 12h, by discrete differentiation of /exp, 

/exp(R, tc + At) - /exp(R, tc - At) 


Respite I R) = 


2At 


The error bars in Eigs. 2C-E are calculated from according to 


APexp(ic I R) 




1/2 


4At2 


Theoretical frequencies are computed using Eq. We find reasonable agreement within 

bars. 


(36) 


(37) 


error 
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3.3 Mean closure time as a mean first-exit time 

The mean closure time {tc{R)) is defined as the mean first-exit time to zero starting from an initial 
radius r(0) = R [30]. We show that it admits a simple analytical expression for the stochastic 
process defined by the Langevin equation § with an absorbing boundary at r = 0 and a reflecting 
boundary at r = i?. 

The solution of the differential equation 


supplemented with the boundary conditions Ti{x = 0) = 0, ^Ti{x = R) = 0 reads 


Ti{x) 


1 / 
1-D [Y 


1 + l/D ) 


(39) 


This yields the mean closure time as a function of initial radius {tc{R)) = Ti{x = R): 

This prediction is compared to experimental data in Fig. 2F. 

In the absence of force fluctuations, Eq. ^ reduces to r{t) = — with the solution r{t)‘^ = 
R‘^ — 2^t. For a given initial radius, the deterministic closure time reads 


tf\R) = 




(41) 


As a consequence the ratio 


tr(E) _ D 
{tc{R)) 7 


(42) 


is always larger than unity: in the presence of fluctuations, the mean closure time {tc{R)) is always 
shorter than the deterministic closure time t^^^{R). 


3.4 Fluctuations of the closure time 

Higher moments can be calculated iteratively [30]. In the case of the second moment we 

solve the differential equation: 


7 d 
X dx 


T2{x)^D 




T 2 (x) =- 2 Ti(x) 


(43) 


with the boundary conditions T 2 (x = 0) =0 and 
second moment {t^{R)) = T 2 {x = R) reads 


-^T 2 {x = i?) = 0, and expression (39). 


The 




5 + 7/d 

(3 + 7/D) (1 + 7/D)2 4^ 


(44) 


The variance simplifies to 


7(i?) = {tliR)) - {RiR)f = 


2D 

3D + 7 


{tciR)f 


(45) 


and the coefficient of variation of the closure time, defined as the 
to the mean value, is a constant 


atSR) _( 2 D y 
YiR)) \3D + jJ 


ratio of the standard deviation 


(46) 


This prediction is compared to experimental data in Fig. 2G. 
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4 Initial mean square deviation 


In order to explain the diffusive behavior of the mean square deviation observed at short time (See 
Fig. 21) we define 

Y(t) = 1 ( 2 )^ ( 47 ) 

and obtain by substitution in (j^, the Langevin equation for Y{t): 

Y = Di+ 7]{t) (48) 

with a drift coefficient Di = The distribution pc,(Y,t) = p{Y,t \ 0,0) obeys the (forward) 

Fokker-Planck equation 




Introducing the scaling variable 


Vi 


and assuming that Po{Y,t) = G{Z), Eq. (1^ becomes 


j^ + Z- + 0 = 2W-[I.,01 


The scaling Ansatz (50) is thus valid in the limit 

2ViDi < ^ 

where the differential equation obeyed by G{Z) simplifies to 

d^G ^ dG ^ „ 

d^+^dZ+^ = 0 


(49) 

(50) 

(51) 

(52) 

(53) 


With the boundary condition lim^^oo G(Z) = 0, the normalized solution of Eq. (53) is G{Z) = 
' - Given the second moment of Z 


+ 00 

{ZV = I i 


Z^ G{Z) dZ = 1 


(54) 


we obtain 


((r(0) - r(t))^) = 2D (F^) = 2D {Z^) t = 2Dt 


(55) 


as expected for simple diffusion. 

Since the scaling variable is typically close to Z ~ 1 according to Eq. ( [^ , condition (52) 
amounts to t y(^)^. During the early stage of closure, we expect that r{t) r(0), and for 
30 < r( 0 ) < 50/im, using the numerical values ( T7p8 ) of 7 and D, the above inequality will be 
respected when t <C 7 h. Eig. 21 confirms that the approximate scaling solution deriving from 
G{Z) is indeed relevant for t < 2 h. 

Note that the scaling (55) does not depend on a specific functional form of the drift Di(r, t), 
and remains valid for the stochastic process defined in Sec. | 2 . 2 | with a non-zero epithelial tension 


cr, although in a range (52) that depends on the drift coefficient. 


When the diffusion coefficient D 2 (r, t) depends on r as in Sec. |2.3[ we checked that the same 
diffusive scaling also holds for short time t and small deviations r( 0 ) — r(t)\ the r-dependence of 
D 2 can be transformed away by an appropriate definition of Y that generalizes (47). However, 


the factor 2D in Eq. (55) is then replaced by a coefficient that depends upon the experimental 
distribution of radii. 
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